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ABSTRACT 


by 


Consider  the  two-factor  crossed  components  of  variance  model  given 


YiJ  =  M  +  Ai  +  Tj  +  eiJ  for  1  =  1»  I;  S  =  1*  J> 


The  random  variables  A^ ,  Tj ,  e^j  are  normal  and  independent  with  means  zero 
and  variances  a£,  o£.  Approximate  confidence  intervals  are  presented  and 

evaluated  for  the  ratios  pa  *  +  °t  +  °e^  and  PT  =  aT^aA  +  aT  +  aP‘ 


2 


1 .  Introduction. 

Consider  the  twu-factor  crossed  Components-of-Variance  model 
given  by 

YU  '  A1  +  Tj  +  £ij  • 

The  are  independent  unobservable  random  variables  and 

A^dKO.o/)  Tj^NtO.o/)  €ij  ~  N(0,°e2) 

H  is  an  unknown  parameter,  and  the  Y. .  are  observable  random 

*  w 

variables. 

The  analysis  of  variance  table  for  this  model  is  given  in  Table  l. 

Table  1. 
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—  ?  2  2 

The  random  variables  Y,  S.j  Sg  .  are  complete  sufficient 

2  ? 

statistics  for  this  problem  and  /e^~  x  (n^for  i = 1 ,2,3  and 
are  mutually  independent. 

An  important  problem  in  applied  statistics  is  to  obtain  confidence 
Intervals  on  various  functions  of  the  variance  components  o^2,  o.j.2  ,  and 

For  some  well-known  functions  exact-size  confidence  intervals  a*e 
available,  but  in  many  situations  exact-size  confidence  intervals  have 
not  as  yet  been  derived. 

For  the  following  important  and  useful  functions  exact-size 
confidence  Intervals  are  known  (exact-size  means  exact  confidence 
coefficient):  a2,  **/<*£,  o_Z/oT2,  4  °A2^’  4  °t2)» 

o//(°£2  +  OpL) ,  o^/{is2  +  o^2).  For  a  discussion  of  these  see  Greybiii  (1976). 

Very  good  approximate  confidence  intervals  are  known  for  o A 
and  Oj2 

For  approximate  confidence  intervals  on  linear  combinations  of 
2  2  2 

°A  ’  °T  *  °e  see  We^c*1  (1950). 


In  applied  work  it  is  often  important  to  determine  the  ratios  of 

2  2  2 

the  Individual  variances  ,  Oj  ,  and  to  the  total  variation 

2  2  2 

( o ^  +  Oj  +  o  ),  1.e.»  to  determine 

*  2 .1  2.  2.  2\ 

PA  °A  /(oA  4  °T  +  °c  } 

.  Z/(  2  .  2  .  2. 

PT  °T  ' '°A  4  °T  4  ae  / 

2,,  2  .  2  .  2, 

P  ■  o  /( ofl  +  oT  +o  ) 
e  e  n  i  c 
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There  Is  no  method  available  for  setting  exact  1-a  confidence 
Intervals  on  pa  ,  py  ,  or  p£.  Large  sample  approximations  have 
been  proposed  by  Osborne  and  Paterson  (1952). 

The  purpose  of  this  article  is  to  obtain  good  approximate  1-a 
upper,  lower,  and  two-sided  confidence  intervals  for  p^  and  Py 
and  evaluate  how  close  the  approximate  confidence  coefficients  arc  to  the 
specified  1-a  . 

The  quantities  p^  and  py  are  useful  In  many  applications. 

Graybill  (1976)  gives  an  example  of  a  two-factor  crossed  components  of 
variance  model  without  interaction.  An  investigator  wants  to  examine 
the  effects  of  drivers  and  automobiles  on  gas  mileage.  One  type  of 
automobile  is  used  and  20  are  selected  at  random  from  one  month's 
production.  Ten  people  are  selected  at  random  from  a  large  city  and 
are  asked  to  drive  each  car  for  two  weeks  and  record  the  number  of  miles 
per  gallon  of  gas  used  for  each  car  for  the  two-week  period.  The  total 
variation  In  the  number  of  miles  per  gallon  of  gas  used  for  this  type  of 
automobile  is  °j\  +  °f  +  °e’  where  is  the  variation  due  to  the  cars, 
oT  is  the  variation  due  to  the  drivers,  and  a?  is  the  variation  due  to 
all  other  noncontrolled  factors.  The  proportion  of  the  variation  in  the 
number  of  miles  per  gallon  due  to  the  cars  is  given  by  ,  and  the 
proportion  of  the  variation  in  the  number  of  miles  per  gallon  due  to  the 


drivers  is  given  by  py. 
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2 •  Approximate  Confidence  Intervals. 

In  this  section  approximate  confidence  intervals  on  pA  and 
Py  will  be  derived. 

First  consider  a  1-a  upper  confidence  interal  on  pA 
denoted  by  *A  <  oA  where  i_A  is  the  lower  confidence  point. 

It  is  simpler  to  find  a  1-a  upper  confidence  interval  for 

a  -  J°A2  „  61  ~  93 
A  I(oTZ  ♦  oc2)  e2  +  "le3 

and  convert  this  to  a  confidence  interval  for  pa  by  using  the  fact 


(  IeA_ 
V  +  f°A, 


(2.1) 


A  function  gL(Y)  of  the  observable  random  vector  must  be 


determined  such  that 


P[9l(Y)  <  eA] 


Is  close  to  a  specified  1-a  for  all  values  of  the  unknown  parameters. 
Since  Y,  ,  S2  ,  S3  is  a  set  of  complete  sufficient  statistics 
for  this  problem,  the  lower  confidence  point  should  be  a  function  of 
these  random  variables.  So  a  function  hL(Y,S^  ,Sg  ,S3  )  must  be 
determined  such  that 


P[hL(^,s12,s22,s32)  <  eA] 


Is  close  to  a  speclfieo  1-a  . 

If  any  constant  k  is  added  to  each  observation  Y..  the 

*  J 

2  2  2 

parameter  oA  is  unchanged  so  it  is  required  that  Iy(Y+k,S^  ) 
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remains  unchanged.  The  constant  k  =  -  Y  is  chosen  and  thus 

hL(T,S12,S22,S32)  =  hL(0,S12,S22,S32)  =  qL(S12,S22,S32),  which  is  a 
2  2  2 

function  of  ,S2  ,  and  S3  only. 

If  each  observation  Y.j  is  multiplied  by  a  non-zero  constant 

c  ,  the  parameter  6^  is  unchanged  so  it  is  required  that 
2  2  2 

q^(S]  ,S2  ,S3  )  remains  unchanged  also.  The  constant  chosen  is 


2  2  2  2 

which  is  a  function  of  /S3  =  and  S2  /S3  =  F2  only. 


The  problem  is  reduced  to  determining  a  function  L(Fj,F2)  such 

that 

P[L(F1,F2)  <  eA] 

will  be  close  to  a  specified  confidence  coefficient  1-a  . 

To  determine  the  function  L(F^,F2)  the  following  five  conditions 
are  imposed. 

1)  The  confidence  interval  must  be  exact  when  e3  =  0. 

2 

When  e3  =  0  it  follows  that  S3  =  0  with  probability  one, 

F-j-^  =  0  with  probability  one,  and  =  0j/e2  •  S°  an  exact  1-a 
lower  confidence  limit  on  is 


F  „  S, 
a:n^  ,n2  c 

and  L(F^,F2)  must  be  equal  to 


In  order  to  satisfy  this  condition. 

2 

2)  When  -*»  it  is  required  that  the  confidence  coefficient 
2  2  2 

be  exact.  When  ■**  and  Oj  ,o£  are  fixed  it  follows  that 

6^  -*®  with  ©2’83  ^xed»  and  e]  dominates  .  So  the  confidence 


Interval  on  should  behave  like  an  exact  1-a  confidence  interval 

2  -1 

on  6j  .  So  L(Fj,F2)  should  be  equal  to  S^  Fa.n  „ 

•  -  1  * 


3)  The  confidence  interval  is  required  to  be  exact  when 
J  -*»  and  I  is  fixed. 

2 

When  J  -*»  it  follows  that  n2  •*»>  n3  -*»  and  hence  S2  -*■  e2 
In  probability  and  S32  -*  e3  in  probability;  also, 

(S22  +  n^S,2)  -  (fl2  +  0,63) 

In  probability. 

Begin  with  an  exact  (1-a)  confidence  interval  for  ©1  given 
by 

S,2F  <  0, 

1  a:n-j  —  1 

and  subtract  S32  from  the  left  side  and  the  "equivalent"  value  e3 


8 


from  the  right  side  to  obtain 


S  2f  -1 
51  a:n 


1*" 


0 


3 


Then  divide  the  left  and  right  sides,  respectively,  by 
2 

+  n^S3  )  and  the  "equivalent"  value  (©2  +  n^)  to  obtain 


s  2f  -1 


1  a:n^ ,«  3  e1  -  e. 


S2J  *  "lS32 


e2  +  nl®3 


Hence,  when  J-*»  the  lower  confidence  point  should  behave  like 

L<FrF2>  ■  -F^-rV 


4)  When  I-+®  and  J  is  fixed  the  confidence  interval  is 
required  to  be  exact. 

For  I-*”  and  J  fixed  it  follows  that  n-j-*®  ,  n3>®  and 

2  2 
S.j  -*■  e1  in  probability,  and  S3  -»  03  in  probability;  also, 

2  2 

S,  -  S3  •+  e-j  -  e3  in  probability. 

Begin  with  a  1-a  exact  confidence  interval  for  e2  given  by 


09  <  So  F  .  „ 

2—2  a:®,n. 


Add  n^S3  to  the  right  side  and  the  "equivalent"  value  n^e3  to  the 
left  side  to  obtain 

0 o  +  n,e,  <  S  2F  .  +  n,  S»2 

2  13—2  a:®,n2  1  3 


Take  the  inverse  of  each  side,  then  multiply  the  right  and  left  sides 

2  ? 

respectively,  by  S-|  -  S3  and. its  "equivalent"  value  -  c3 


9 


to  obtain 


S  2  -  S  2 


61  "  e3 


S  2F  +  n  S  + 

*2  V«,n2  n1S3  2  1  3 


Hence,  when  I-*~,L(F^  .Fg)  is  required  to  behave  like 


F,  -  1 


L(FTF2)  =  f2f  +  n* 

2  a:«,n,  1 


5)  If  an  a-level  test  of  H0:e^  =  g3  vs.  >  e3  is 

accepted  the  confidence  interval  should  include  "zero"  and  L(F^,F2) 
should  be  increasing  in  for  fixed  values  of  F2  . 

H  is  accepted  iff  F,  <  F  ,  so  the  confidence  limit 

u  i  —  a . n>j  ,n3 

should  be 

L(F,,F2).0  when  F,  <  F^  ^, 

L(FrF2)  monotonic  increasing  in  F-j  for  fixed  values  of  F2 
when  F,  >  F 

I  a:nj,n3  . 

The  form  of  the  function  that  could  be  used  for  L(F^,F2)  is 

a„  +  3-iFt 

L(FrF2)  "  b  +  b,F? 

O  1  2 

but  a  more  general  function  will  be  examined,  i.e., 


L(FrF2) 


.  ao  +  alFl  +  a2fl 


bQ  +  b,F2  +  b2F2‘ 


where  a.,  and  b^  which  are  functions  of  n^,n2,n3,n  are  determined 
by  requiring  L(F1,F2)  to  satisfy  conditions  1),  2),  3),  4),  5) 
Substituting  condition  1)  gives 


-  F 


ct .  n  i  s  n  ^ 


(2.2) 
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Substituting  condition  2)  gives 


and  hence 


From  condition 


b,  =  F  F 
1  a:n^ ,»  a:n^ ,n2 


3)  Lim  a  =  -1  ,  Lim  b 
J-xo  0  J-x»  c 


Lim  a, 
J-XO  1 


(2.3) 

(2.4) 

F 

a:n1 


Lim  b,  =1,  Lim  a2  =  0  *  Lim  b2  =  0 

J-K»  J*X»  J-X» 


which  are  satisfied  for  those  values  of  a-j ,  b-j  from  equations 
(2.3)  and  (2.4). 

From  condition  4)  Lim  a  =  -  I ,  Lim  b  =  n,  ,  Lim  a,  =  1  » 

I-**  0  I-x»  0  1  I  -Xo  1 


Lim  b,  =  F  „  ,  Lim  a0  =  0  ,  Lim  b,  =  0  . 

I-**.  1  a:“»  n2  I-xo  2  I-~>  2 

Hence  the  limits  of  a^  and  b^  are  consistent  for  condtions  3) 
and  4). 

From  condition  5)  L(F1>F„)  =0  if  F,  <  F  _  which  implies 

'12  1  —  o:n, ,n_ 

that 


an  +  a-iF,  +  a5F,“^  =0  if  F,  <  F  . 

o  li  2  1  l  —  a:n.j  .n^ 


Now  substitute  F,  =  F  into 

I  Q • Di 


-l 


a„  +  a,F  _  +  a,F  „  =  0  , 

0  1  C(  •  0|  2  Cl •  Di  j 

let  a-j  be  the  value  given  in  equation  (2.3)  and  solve  for  a2  to 
obtain 


a2  *  <-»„ 


-  F 


F 

a.n^  ,qo  ct •  n i  > 


)F 


a • n i > 


.  j 
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Use  the  simplest  values  for  aQ  and  bQ  ,  namely,  aQ  =  -1 

and  b  ®  n,  ,  and  obtain 
o.l 

-1  +  FalL-  F1  +  (1-FalJ,  -F«:ni,n-,F«:ni  n,  F1 

L(F  F?)  -  - I - T - L 

'  n,  +  F  "I  F  _  F0  +  b„F0” 


-1 


1  a:n^>®  a:n^ ,n^  2  2  2 


If  F1  >  F  n  n 

1  a  H| ,n^ 


(2.5) 


L(F,,FJ  =  0  If  F,  <  F  .n  n 
1  Z  I  —  a:n-j  ,n. 


To  determine  the  value  of  b^  we  impose  condition  5). 

2  2 

Condition  6;,  When  HQ:oy  =0  vs.  H^Cy  >0  is  accepted 

the  confidence  interval  is  required  to  be  exact. 

2  2 

To  examine  condition  7)  we  get  HQ:oy  =0  vs.  H^:oy  >0  <“> 

*;02  =  e3  VS-  Hi:e2  >  e3  w  Ho:9A  =  (ere3)/I03  vs- 

H1:0A  <  (°r03)/l03  *  /e  v 

If  HQ  is  true  then  eA  =  Jo^/Io^2  =  f  -  1  j 

So  if  Hq  is  accepted  at  a  level  we  want  L(F^,F2)  to  be  an 


H 


exact  1-a  lower  confidence  point  for  eA 

H  is  accepted  iff  F,  <  F  „  _  .  An  exact  l-o  lower 

o  r  2  —  a^in^ 

confidence  point  for  -  1/1  is 


s  2  F  -1 

1  a:nl,n3  1_  91  '  93 

IS,2  1  ~  1  °3 


So  when  F0  <  F  . 


2  —  a:n2’n3  we  want 


s  2f  -1 

r~., 


L(F 


F  x  -  1  a:rVf13  1 

r2)  ‘  IS2  '_I  * 


This  is  not  consistent  with  the  fir*,  live  coik.  :  t  ioo  .  so  we  only  require 
the  denominator  of  (2.5)  to  be  I  when  condition  <>)  is  satisfied. 


12 


So  set 


and  then 


nl  +  ^a:n^ ,®^o:n^ .ng^ccng.n^  +  ^a-.ng.n^  “  * 


^2  ^  ^a:n^  ,®^a: n^  ,n2^a: ng .n^^a: n2 .n^ 


and  this  satisfies  Lim  b2  =  0 

J-*vO 

The  final  result  is 


-1  +  F  F,  +  (1-F  “If.  )F  f:1 

L(F  f2)  = - TT - i - A - - rr 

»  c  n+F*  F  f  +(1  -F  F  F  \  r  F* 

1  a.*n|,®  a:n^  ,r\^  2  '  a:n^,®  an-j  ^  a:n2»ng  a:n2»n^  2 


if 


's/a/if 

■  J  a.n1 ,n3 
and 

S  2/S  2  >  F 
2  3  a. ngiHj 


L(FrF2)  = 


-i  +  f  :I  f,  +  (i-F  ;!  f  )f  n  f: 

a:n^,®  l  '  a:n^,«  ain^.n^  a:n-j,n^  1 


-1 


(2.6) 


/ 


if 


S,2/^2  >  F 

1  3  a.n, ,n~ 

and 
2  2 

SJ VS,  <  F„.n  _ 

2  3 


L(F,,F2)  =  0  if  S,ViF.:Vn, 

I  V 
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Instead  of  using  L(F-j  .F^)  in  equation  (2.6),  the  simpler  function 
in  equation  (2.7)  will  be  used;  which  is  obtained  from  equation  (2.5) 
by  setting  b2  =  0  .  Note  that  Lim  b£  in  equation  (2.6)  satisfies 

I-Ko 

the  requirement  in  condition  5).  So  the  lower  confidence  point  for  is 


-i  +  r  7*  f,  +  (1-f  _1  f  )f  f; 

a:n^,«  1  '  a:n^,»  a.-n^.n.,'  a:n^,n3  1 


la(fi»f2) 


r-1 


1  a.rij  ct •  n i  > n 2  2 


if  S  2/S  2  >  F 

1  3  a:n^,n^ 

2,e  2 


LA(F,.F2)-0  If  h/S3lf„:nrn 


(2.7) 


Note  that  the  confidence  interval  in  equation  (2.7)  covers  the 
maximum  likelihood  estimator  of  0^  ,  say  eA  .  To  demonstrate  this 

A 

the  M.L.  estimator  of  oA  ,  denoted  by  9A  is 


eA  =  (Fi"1 )/(F2+nl )  if  Fi  >  1 
eA  =  0  if  F1  <  1 


> 


and  the  confidence  interval  is 


-1  +  F  "l  F,  +  (1-F  F  )F  F,_1 

a:n^,«  I  a:n^  ,<«  ain^.n^  ain-j.n^  1 


nl  F«:n1.-F.:n1.n2F2 


—  eA  <- 


if  F,  >  F  n 
1  o  •  n^  fii^ 


0  <  6.  <“ 
-  A 


if  Ft  <  F  .n  n 
1  —  a :  n^  .n^ 
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Now 


•1  ♦  F„:I  .F,  ♦  O-F  '1  J 


)F...  „  p;'  <  F,  -  1  , 


ainp*  1  v  a:^/0  ain^nj'  ain^rij  1  1 


and 


Mn,,"*  ’  • 


Also 


(1-F  F  )  <  0. 

'  a:n^  ,oo  a:n^  tr\^' 


nl  +  ^a:n^  ,«o^o:n^  ,02^2  ’  F2  +  nl 


because 

Fa:nj ,«Fa:n.| ,n2  >  ^ 

*  A 

So  *"A'^1  *F 2^  <  eA  <  “  anc*  0a  ^  in  the  l"a  confidence  Interval. 

Next  consider  a  lower  l-o  confidence  interval  for  eA  .  A  lower 
1-a  confidence  interval  for  oA>  given  by  0  £  oA  5  UA^rrF2^  can  be 
readily  obtained  by  replacing  a  by  1-a  in  equation  (2.7).  The 
upper  1-a  confidence  point  of  the  lower  confidence  interval  is 
given  by 


■1  +  F 


-1 


-1 


ua(fi»f2) 


l-a:n^,®Fl  +  ^  Fl-a:n^  ,«°Fl-a:ni  ,n^Fl-a:n^  ,n.jFl 


-1 


Hi  +  F1 F,  F„ 

I  l-a:ni ,®  l-ain^ ,n2  2 


VF1’F2>  =  0 


If  F1  >  F1  n 

I  l-a.n^  * 

If  F1  1  F1  n 

I  I -a : n^ .n^ 


(2.8) 


The  procedure  for  determining  this  lower  confidence  point  is 
similar  to  that  described  for  the  upper  confidence  point. 

From  the  formulas  for  upper  and  lower  confidence  intervals  in 
equation  (2.7)  and  (2.8),  the  two-sided  confidence  interval  on  r. 
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can  be  readily  obtained.  The  1-a  two-sided  confidence  interval  for 


is  given  by  LJ  <_  <_  ,  where 

.-1 


LX 


~}  *  Fa/2:nr-F1  +  (1"Fa/2:n1  ,-Fa/2:n1,n3)Fa/2:n1  .n^i* 

™  «  "  a  C  tr  r 


n1  +  Fa/2:n1rFa/2:n1,n2F2 


(2.9) 


UX 


-1  +  Fl-a/2:nr»Fl  +  (1"F1-l/2:n1,-F1-a/2:n1,r.3)Fl-a/2:n1,n3Fl1 

nl  +  F1-a/2:n1.«F1-a/2:n1,n2F2 


This  will  be  proved  in  Chapter  3. 

If  either  limit  is  negative  it  Is  replaced  by  i.ero. 

As  mentioned  in  Chapter  1,  due  to  the  symmetry  in  the  model.  If 

a  confidence  interval  fcr  is  obtained,  the  confidence  interval 

for  Oj  can  be  easily  obtained.  The  formulas  for  the  confidence 

7  7  7 

Intervals  for  eT  are  given  below  where  Oj  =  Io^  /(o^  +o£  )J.  The 
l-o  lower  confidence  point  for  0^  is 


-1  +  F  F,  +  (1-F'!  F 

o:n0,»  2  ' 


lT(FrF2) 


)F... 


--1 


n,  +  F  ."J  F  F. 

c  a:n2 ,«  a:n2 ,n^  1 


o:n2,“  a:n2,n3;  a:n2,n3  2 


(2.10) 


The  1-a  upper  confidence  point  for  eT  is 


UT(FrF2)  =  - 


*  FNa;„„-F2  *  ,  F 


)f, _  -  f: 


2 1  -a.r^n^  1  *ci inp 2 


n,  +  F/  F, f7 
2  l-a:n2,~  l-a:n2,n^  1 


-  (2.1D 


The  1-a  two-sided  confidence  points  for  are 

-1  +  F  ,1  F,  +  ( 1  — F  .1  F  ,,,  _  )F  /9.„  „ 

a/2:n2*®  2  a/Z.n^i®  a/2:n2»n2  a/2:n2,n2  2 

.1  "*  *  f«/2:n2,.f./2:n2,n,Fl  (2 

^  *  Fl -a/2: n. .—F2  ^  F1  -a/ 2 : ,  f. ^ ^  2 

LT  "  T} 

n2  +  F1-a/2:n2,»Fl-a/2:n2,n1Fl 

If  any. limit  is  negative  it  is  replaced  with  zero. 

From  the  confidence  intervals  obtained  for  and  Oj  »  and  by 
equation  (2.1)  we  get  the  corresponding  confidence  intervals  for  pa 
and  pj  . 

The  1-a  lower  confidence  interval  for  p^  is 


(o  +  aA/ 

where  LA  =  L^F^)  wh1ch  is  defined  in  equation  (2.7). 
The  1-a  upper  confidence  intervals  for  pa  is 


(J3*  \ 


<  Pa  <  * 


where  UA  =  ua^ F1 »F2 ^  which  is  defined  1n  equation  (2.8). 
The  1-a  two-sided  confidence  interval  for  p.  is 


\  l  H 

-  pa  -  \j  ♦  il”a 


where  LJ  and  are  defined  in  equation  (2.9). 
The  1-a  lower  confidence  interval  for  py  is 


V 
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where  Lj  =  Lj(F^ ,Fg)  which  Is  defined  In  equation  (2.10). 
The  1-a  upper  confidence  Interval  for  pT  is 


where  Uy  =  UyfF^.Fg)  which  Is  defined  In  equation  (2.11). 

The  1-u  two-sided  confidence  interval  for  ov  is 


i  luuiicc  mtci  va  i  iui  My 

(az).  I'h-) 

\J  +  lUf  J  -  PT  1  \J  +  ILj ) 


Where  and  Ly  are  the  confidence  points  given  in  the  equation 

(2.12). 
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3.  Evaluation  of  the  Methods 


In  the  preceding  section  confidence  point(s)  were  determined 
for  upper,  lower  and  two-sided  confidence  intervals  which  have  exact  an 
confidence  coefficient  1-a  for  specified  restrictions,  and  for  general 
conditions  the  confidence  coefficient  should  be  close  to  1-a  .  In 

this  section  methods  will  be  derived  to  determine  how  close 
the  approximate  confidence  coefficients  are  to  the  specified  ones. 

First  consider  the  upper  confidence  interval 

PCL(F1  ,F2)  <  oA  <  »]  =  PL(eA)  (3.1) 

From  equation  (2.7) 


-1+F  F  +(1-F  F  IF  F~^ 

Sam-.,”  1  ain-.,”  ain^n,'  a:n,,n.  1 

- ! - U. - L_J - L_J -  if  F.  >  F  .n  n 

n  x  c  *  p  p  ci  •  Hi  »n, 

1  a:nj ,®  a:n^ .n^  2 


If  F,  <  F  n 

1  ot • n i 9 

(3.2) 


From  Section  1 


F1  "  Si2/S32  for  1  "  1*2  i  als0  n^si2/ei  =  ui 


for  1  =  1,2,3  are  Independent  chi-square  random  variables  with  n. 
degrees  of  freedom  respectively. 

Consider  the  following  events 


A  =  (L(Fj,F2)  <  eA;  F1  >0;  F2  >  0} 


"  <F1  >  Fa:n,,n3'  F2  >  "•  L<FrF2>  ~  0! 


.1 
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»  <ft  <  Fa:ni,n3»  F2>0’  ^^>0}  '  (3.3) 


So  B^n  Bg  B  0  and  U  B2  *  n  where 


«■  |f,.F2.UFrF2))  :  F,  >0i  F2  >  0i  L(F,.F2)  >  oj  . 

Clearly  PL(eA)  s  P(A).  We  will  write  events  and  B2  In  (3.3)  as 

B,  ■  <F,  >  F  n  }  ;  B9  =  tF,  <  F  „  }  (3.4) 

1  1  a:n^,n^  2  1  —  ain^.n^  '  ' 


From  this  we  obtain 


PL(eA)  *  P(A|B1)P(B])  +  P(A|B2)P(B2) 

W  -  PCLCF, ,F2)  «  «a|F,  .  Fo;ni>n33F[F,  >  Fo;„]in3] 


♦PWF,^)..*^  <  F.^MF,  <Fo:nj>nj3 

”1  c  x/i  r  "1  c  \r  r”l 


PL<eA>  ■  p 


1+F 


a:n.j  ,»»F1+^  Fa:n.j  ,*Fa:n^  ,n3^Fa:n.|  ,n3Fl 

n-.  +  F  ”  F  F0 

1  a:nr-  a:nrn2  2 


ieA 


F1>Fa:n1 ,n3 


P£F1  >  Fa:nrn33  ♦  P[0  £  eA|F,  <  F^^PCF,  <  F^^] 


But  eA  Is  always  greater  or  equal  to  zero  so 
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PD.(F,,F2)  <  eA|  F,  iF„;niinj].l 


W 


i  ,  p  “1  C  x/  1  p  **  1  r  \r  r* 

ro;n1>«  1  v  a:n^,*>ra:n^  ,rig  1 


nl  +  Fa:n^ ,«Fa:nj ,1x^2 


-  0A  F1>Fa:n1,n3 


<3-5> 


The  first  term  in  (3.5)  gives 


"-1+F  ^  F  +(1-F”^  F  )F  F”^ 

o:n, ,«  1  '  a:n-|t<»  a:n,  ,n,  a:n,  ,n,  1 

p  =  p  - ! - - - L_3 - <6  F  >F 

1  n,  +  F  f;  F  F?  A  1  a,nrn3 

1  a.n^  *°°  d . n i  > n ^  2 

Fa:n,  ,®+Fl+'Fa:n, ,<*>  Fa:n,  ,n,^Fa:n,  ,n,Fl 

-  P  - ! - ! - L-J - \ — i - <  0  F  >f  n 

n,F  +  F  .  F,  A  1  °-Vn3 

L  1  a.n^ ,«  a:n^  jng  2 


pi  - p 


F  S,2$,2+sA(F  -F  )F  „  S  4 

ain-j,®  l  o  I  a:n.|,°°  a:n^,n3'  a:n.j,n3  3 


niF  •„  S,  S,  +  F  .n  „  S/S  ‘ 
1  q<  n-j ) °°  3  1  ci«  H|  *^2  ^  • 


i8AlfrF.;n,.n3 


2P  2  .  P  4 


Pi  *  p[-F  S/S.  +  S,  +  (F  -F  )F  S  < 

1  L  a:n-j,»  1  o  1  '  a:nj,®  a:n^  .n^'  a:n^  ,n3  3  — 


n,F  S2S,20A  +  F  .n  n  s/s/o.jF,  >  F  ,n  1 
1  a:n^,«  1  3  A  ain^n^  Z  1  AM  a.n^.n^-* 


Pi .  ,  (-AA  m  ,  (F  ,  )F  feVW 

\  el2  /  nl  “-"l-  “:Vn3  °:V"3\  <,//  „32- 
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nlSl  Yn3S3  |  9103 


F.:n,  ,J1+nl°A>  H77I  03-/  njn^  +  ^  .n^A  [ TJ 


p  >  p 

V  J  nln2  ^  a:n^ 


Ul2el2  t^e  2 

P  3  p  -i  +  ( F  -F  )F  *2  <  c 

1  2  '  c*:n, a.*n,  .n,  a:n,  ,n,  _  2  —  a:n,,  (1+n.e.) 

n-j  1  13  1  3  1I3  1  «'  1  A 


®1®3 

x  Mo  F 


ele2 

8«  M0  ~~  F,  >  F 


r 3  n]n3  T  ra:nrn2°A  °lu2  n1n2  rl  Vn^ 


where  =  n^sZ/e^  for  1  =  1,2,3. 

o 

Divide  both  sides  of  the  inequality  by  e3  and  let  p..  =  0^/63 
for  1  ■  1,2  to  obtain 


P]  -  P 


11  2  2 
U1  P1 


72”  +  (Fa:n1..'Fa:n1,n3)Fa:n1,n372  -  Fa:nr»(1+nleA) 
1  3 


U1U3 

*  nn  Pi  +  F  • 


MJ  F,  >  F 


n^n3  H  a:n^,n2  A  n^n2  1  2  1  a:n^,n. 


Divide  both  sides  of  the  inequality  by  U-j  to  obtain 


V  ^  V  ^,r.-r.:V,3>  ^ 


F°:nrn3  < 

P  2  U  ^  a !  n-j ,  oo(  1 +  H-J  0^) 


«  U3Pl  4.  r  0  plp2  .a  p  p 

a.n^fn2  ^^1^2  ^  *  ci  •  n  ^ 


2  F  2 

Pi  #  o:ni,n,  U,fc 

^2  U1  +  (Fa:n1,«-Fa:n1,n3)  — 2~  '  F.:ni,-<1+VA> 


plU3\  nln2 


n1n3  /  0A°1P2  Fa:nrn2  1  T  *  ^r»3 


p  "?P1 _ u  +  (Fa:nr<n"Fa:n1»n3)Fa:n1>n3 

W77;.:n1,n2  1  n3plp2eA  Fa:nrn2  '  'U1 


(1+nl°A)Fa:nr- 
nlp2eA  fa:n]tr\2 


U,  i  U,  F,  >  F  in  n 
<5  l  1  a.n-j  >n^ 


Let  ci  =  (n2pl)/(nl8Ap2  Fa:n1,n2) 

C2~  *Fa:n1,-'Fa:n1,n3lFa:n1,ni  (n3plp20AFa:n1  ,n2^  (3*6^ 

C,  =  (1+11,6.  )F  , 

J  I  A  a:n^,y  1^2  A  a:n^,n2) 

Note  C-j  >  0  ,  C2  <  0  and  C3  >  0.  Then 

u/ 

VP[C  )UltC2-^-C3U3<U2|F1,Fa;n)tn3] 
[c,u,(u,*u3)  C2u32  (U,*U3)2  C3u3 

'  [W  777  “7“  '  TVU37  1 W  i 


^*2  F1  >  Fa:n^  ,n3 


1 


h 


Then 


d 

11 
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,  fa:n.  ,n,+  n2°l 
(1-Y)'1  >  -  1  3 


a.n-j 


F  _  +  n,p, 
a:n^ ,n3  21 


(3.8) 


n2pl 

Then  (3.7)  can  be  written  as 
,PL(eA)  =[(C1Y  +  C20-Y)2Y_1  -  C3(1-Y))x  <  U2  ; 

Y  >  Fa:n^  .n^/^ct.-n^  ,n3+n2pl 

+  P|Y  <  F  n  /(F  „  +n0p,)l 

L  —  o:npn^  a:n^,n3  2  1  J 

=  P  [((Cj+Cg+C^Y+CgY"1  -  (2C2+C3))X  <  U2;  Y  >  a]  +  P[Y  <  A] 


where 


=  P[g(Y)X  <  U2;  Y  >  A]  +  P[Y  <  A] 

9(Y)  =  (C]+C2+C3)Y  +  C2Y_1  -  (2C2+C3) 

X  e  F  /(F  +  n0p, )  . 

o. n-j  )D^  ct  •  n -j  jp j  2  1 


(3.9) 


With  probability  one  Y  e  (0,1]. 

Consider  the  function 

(0,+C/Cjy2  -  (2C?+C  )y  +  C? 

9(y)  -  - 1  2  -3- — - ^ - £  0  <  y  <  1 


9(y)  =  0  iff  (C1+C2+C3)y^  -  (2C2+C3)y  +  C2  =  0; 


Clearly 


hence,  g(y)  has  only  two  zeros  for  ye  R,  and  they  are 


.  (X2*C3)  *  [C32-4C,C2]1/2 

.1  ■  2(C^+Cp 

.  '  tC32-4C,C2]1/2 

°2  "  2(C]+C2+C3) 

Also  Urn  g(y)  -  -®  and  Limg(y)  =  C,  >  0.  Since  g(y)  is  continuous 
y-»0+  y-+l 

in  (0,1]  then  at  least  one  of  the  roots  (6^,62)  Is  in  (0,1).  But 
% 

^jM  =  (Ci+C2+C3)  .  c^'2  =  C1+C3+C2(l-y'2)  ,  0  <  y  <  1 

But  for  y  e  (0,1),  ^1  -  *4^  <  0  .  Also,  (C^+C^)  >  0  and  <  0  , 
hence,  if  ye  (0,1)  then 

Mli  >  o 
dy  u  • 

So  g(y)  is  monotonic  increasing  in  (0,1];  also  since 

Lim  g(y)  =  -  •  and  Lim  g(y)  =  C,  >  0  ,  only  one  of  the  two  roots 
y-K)+  y*l 

can  be  in  (0,1].  This  gives  the  following  two  possible  forms  for  g(y). 


Figure  3.2 


The  number  6  is  the  root  in  (0,1);  also  Ae(0,l)  so  if  g ( x )  <  0 
this  implies  that  A  <  6  .  Next  we  show  that  g(A)  <  0  for  some  value 
of  A  €  (0,1). 


g(A)  =  (c1+c2+c3)x  +  c2a''  -  (2c2+c3) 


(Cl+C2+C3)Fa:nrn 
^Fa:n1,n3  +  n2pl^ 


l’n3  A  ^  a:nl  ,n3 


'2  F„. 


a:Vn3 


•  (2C2+C3) 


^Cl+C2+C3)Fa:n1  .n/^a.-n,  ^(F^^+ry, ) 
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Replace  C^,  C^»  and  by  their  values  in  (3.6)  and  obtain 

•  P2  2  2C  P  2  2-2 

n2pla:n,,n,  n2  P1  Vn. ,»  a:n,  ,n,  "  n2  P1  a:n.,n, 

N  *  — - pJ — ^ - +  1  1  J 

nTpo^A ^ 


ip20A  a:nrn2 


n0p , p08.F ~ „ 
3  1  b  A  ct  •  n  j )  n* 


n9  (®Ani  +  ** )PiF  .n  F 
«  A  1  1  d •  n j  j®  ot  •  fi| 

rilpIO  # 

3  2  A  a.n^ 


['^Aa:^ ,n3‘n22pl2Fo:n1 ,n3+Fa:n1 ,*n22pl 2-n22pl 0Anl +1 )Fa:n1 ,-JFa:n1 ,n, 

n3plp2eAFa:nrn2 


Fa:nr«Fa:nrn3n22pl2(1"0Anl+1) 

n3plp20AFa:nlfn2 


■0Anln22pl2Fa:n1,®Fa:n1,n3 
n3plp2eA  Fa:nrn2 


<  0 


So  the  numerator  is  always  negative  which  implies  that  g(x)  will 
always  be  negative;  and  hence,  A  <  6  . 

Summarizing  these  results  give 

a)  g(y)  i  0  iff  y  ^  6  for  y  e  (0,1) 

b)  0  <  A  <  5  <  1  ' ’  (3.10) 

c)  g(y)  <0  iff  y  c  (A, 6) 

Using  these  facts  to  evaluate 

P[g(Y)X  <.  U2,  Y  >  a]  in  (3.9) 
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we  obtain 

*  P£g(Y)X  <  U2;  Y  >  A|A  <  Y  <  6]  P[a  <  Y  <  6] 

+  P[g(Y)X  <  U2;  Y  >  X | Y  >  6]  P[Y  >  5] 

-  1  •  p[x  <  Y  <  5]  +  [g(Y)X  <  U2;  Y  >  Ai  Y  >  5] 

«  P[A  <  Y  <  6]  +  P[g(Y)X  <  U2;  Y  >  6] 

■  Fy(6)  -  Fy(A)  +  P  [g(Y)X  <  U2;  Y  >  6] 

where  ' 

Fy(y)  =  P[Y  <  y] 

From  the  original  probability  in  (3.9)  we  get 

PL(eA)  =  Fy(6)  -  Fy (A)  +  P[g(Y)X  <  U2;  Y  >  6]  +  Fy  (a) 

“  Fy(6)  +  P[g(Y)X  <  U2;  Y  >  6]  (3.11) 

Also 

P[g(Y)X  <  U2;  Y  >  5]  =  P[g(Y)X  +  X  <  U2  +  X;  Y  >  6] 

■  P[g(Y)  4  1  <  --2-  j  X  j  Y  >  6] 

c  p[[g(Y)  +  i]"1  '  Y  >  6]  (3.12) 

Let  Z  =  X/(X4U2)  ;  It  follows  that  Z  is  distributed  as  Beta  with 

/nl+r>3  n2  \ 

parameters  \ — 2 -  '  ~T>  ’  Also,  let  ni  be  an  even  integer  so  that 

n.j/2  =  for  i  =  1,2,3  and  hence  the  are  positive  integers. 

If  X.|  ,  X2,  X^  ’re  independent  chi-square  random 
variables  with  ,  m2,  degrees  of  freedom,  respectively,  and 
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'1  +  X2 


V2  -  X,  ♦  X2  ♦  X, 


X,  *x2 


then  ,  Y?  are  mutually  Independent  Beta  random  variables  with 
parameters  (m^/2 ,  m^/Z)  and  ((m1+m2)/2,m3/2^  respectively. 
From  the  result  in  (3.12)  we  obtain  (let  v  denote  [g(y)  +1] 

pl(0a)  -  p[[g(YHr 1  >  Z;  Y  >  fi]  +  F y(6) 


6  0 


/'/*  B~  (m,  ,m,)  m,-l  m,-l  nu+m,“l  m..-l 

/  »  <'■*)  2  <’-2>  “  dy 


1  >V  -1 

f'r  B  (m,,m-)  m,-l  m,-1  m  +m,-l 

■  fy({)  *J  J  »’  <>->>  2 

6  0 


ro2’1  /m„-A  r  rn9-r-l 
Z  l  %  )(Dr  (-2)  2  dz  dy 
r=0  x  r  ' 


m  1  1  rn.-r-l 

X  B  (ml’m3)r(m2)(_1) 

FY(6)  +  ZQ  B  ( m'1  +m3  ,m2 )  r  ( r + 1 )  r  [  m^rT" 


[g(y)+i] 


■1 


n  law/ 

/y1"1'1  0-y)"3"'/  2  < VWr-?)  dz  dy 
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/I  _  i  _  ,  -(m.+nu+nu-r-l) 

y  1  (1-y)  3  j-taiilill  '  2  3  dy 
*  '  3‘  (m^+m^+m.j-r-l)  * 


_  i  m,-r-l 

2  (-1)  2 


^  .  .  (nij+nyiiyl)! 

+  r5)  - 1  i  I  (m3-  i )  i  (m2-  r- 1 J  i  r  { ( +n»2+m3-  r- 1 J 


J  m'’1  (l-y)1"3  1 


r  p  o- 

J  (m1 


«  [g(y)+i] 


♦iT^+m^-r-l ) 


dy 


(3.13) 


From  the  expression  In  (3.13)  w<?  will  compute  pL(0^)  • 
Next  consider  the  lower  confidence  Interval 

P[*A  <  u(f,.f2)]  ■  Pg(eA) 

From  equation  (2.8)  It  follows  that 

-  f:1  F,  +  (1-F-1 
l-a:ni ,«»  1  v 

U(F, 


(3.14) 


'  +  ^l-a:nj  f“^l  +  ^  Fl-a:nj  ,»Fl-a:n^  ,n^Fl-a:rij  .n^t 

W  * - TTFr - : - 7~ - 

1  l-a:n-j,«>  l-a:nj,n2  2 


U(FrF2)  =  0 


if  S^/S,2  >  F,  ,n  _ 

1  3  l-aOj, 

If  S,2/S  2  <  F,  n 
1  3  1  ”  ou  n  i  y 


(3.15) 


Substitute  equation  (3.15)  into  (3.14)  to  obtain 

W  e  p[Cg(Y)  + 1]_1  -  2>  Y  i  6]  (3.16) 


where  Y  ~  Beta  (m^.m^)  ,  and  Z  ~  Beta  (m^+m^,m2)  ;  Y  and  Z 
are  independent. 
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9 (y)  •  (C1+C2+C3)y+  C2y_1  -  (C3+2C2)  (3.17) 


6  Is  the  largest  root  of  g(y). 

The  C's  in  g(y)  are  obtained  by  changing  the  a's  to 
1-a's  in  the  equation  (3.6).  So  we  have 

n,p. 


r  _ _ £M1 

1  .  n,e.p,F. 


rAp2rl-o:nrn2 


^l-a:n^ ,«  Fl-a:nj  ,n3^l-a:n^  ,n3 
n3p1p26AF1-o:n1,n2 

*1+nleA)Fl-a:nr» 

3  "lp2eA|:l-a:n1,n2 


(3.18) 


Note  Cj  >  0  ,  C3  >  0;  and  C2  >  0  if  n^  >  2,  n2  >_  2,  and 
a  <_  .10. 

Consider  the  function  g(y)  in  equation  (3.17)  .  Clearly 


g(y)  *  0  Iff  (C1+C2+C3)y  -  (2C2+C3 )y  +  C2  *  0; 

hence,  g(y)  has  only  two  zeros  for  y  c  R  and  they  are 

(2C2+C3)  +  [C32-4C1C2]1/Z 

61  *  2(c1  +  c2  +  C3) 


(2C2+C3)  -  [C,^4C,C2] 

*2 " "  Ttf;  7  VY1 


.  -.1/2 


Also  Llin  g(y) 
y-*0+ 


and  Llm  g(y)  ■  C,  >  0 

y*i  1 
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Since  g(y)  is  continuous  in  (0,1],  then  either  both  of  the  two  roots 
of  g(y)  are  in  (0,1]  ,  or  neither  of  the  roots  are  in  the  interval 
(0,1]  ,  or  the  two  roots  are  complex  and  not  real.  But  A  c  (0,1]  so 
that  g(A)  <  0  .  For 

Fl-a:nrn3 


1-a:nj ,n3  +  p1n2 


(3.19) 


9(A) 


2  2 

_0AP1  nln2  Fl-a:n1,-Fl-a:n1,n2 
n3p1p2eA  Fl-a:nrn2 


<  0 


This  follows  from  the  proof  of  (b)  in  Equation  (3.10). 

From  the  above  It  Is  clear  that  the  two  roots  of  g(y)  are  in 
(0,1],  which  results  in  the  following  form  for  g(y). 


From  equation  (3.16)  we  evaluate  Py(e^)  as  Follows 

m  i  nu-r-1 

V'  (_i)  2  (m1+m2+m3-l ) ! 

Pl/eA^  *  (m3-l ) !  (rn^-1 ) !  (r^-r-l )  IrT (rn^+rn^+nij-r-TT 


i 
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/'nu-l  m,-1 

y  1  0-y)  3 

6 

The  two-sided  confidence  Interval  Is 

PU^.Fg)  <  eA  <  U1(F1  .f2)3  -  P2(eA)  (3.21) 

where  L^F^Fg)  Is  given  by  (3.2)  with  a  replaced  by  a/2  and 
U] (Fi ,F’2)  given  by  (3.15)  with  l-o  replaced  by  l-a/2  . 

To  show  (3.21)  let  A  and  B  be  the  following  events 

A  -  {L1(f1,f2)  <  eA)  B  -  {eA  <  U1(f1,f2)} 

Then 

P2(«A)  *  P(Af\B)  (3.22) 

From  (3.22)  we  obtain 

'  WO  B)  -  1  -  P(A7TB)  -  1  -  P(AUB) 

:  1  ~  [P(K)  +  P(B)  -  P(Afl  B)] 

■  1  -  P(*)  -  P(B)  +  P(ff  O  ff)  (3.23) 

Also 

P<*)  -  1  -  PL(eA) ,  P(B)  -  1  -  Pu(eA) 


[g(y)+i] 


(m^+m2+m3-r-t) 


+  1 


dy 
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P(*n  B)  «  p(eA  <  L(F1  ,F2);  U(F1  ,F2)  i  0A)  =  0 

\ 

If  L(F^  ,F2)  probability  one.  Next  we  will  show 

that  this  is  indeed  the  case  for  n-j  >  ?  ,  n2  >.  2  and  a  <_  .10  .  Vie 


must  prove 


+  Fa!ni,-ri  +  (1-Fa!nr-Fo:n1,n3)  Fa:nrn/l 


r*i  *  F  F  F0 

1  a.n*|  i®®  oi •  ^2  2 


^  +  Fl-a:n^  ,°-Fl  +  ^  Fl-a:n^  ,*Fl-a:n.j  ,n3^Fl-a:n-j  ,n3Fl 


nl  +  Fl-o:r.r»Fl-a:nrn2F2 


which  we  write  as 


■1  +  A]F1  +  BjF’1  -1  +  Aj,Fj  +  B^-1 

nl  +  Cl^2  “  nl  +  C2F2 


(3.24) 


This  implies 


U(F1,F2)  <  L(FrF2) 


To  prove  equation  (3.24)  we  use  Lemma  l,  2,  and  3.  (See  Ghosh  (1973)). 
Lenina  l .  F _  „  Is  decreasing  In  n,  if  n,  >  2  ,  n0  >_  2, 

-  o.n^ng  c  i  ^ 

a  <  .10  • 


Lemma  2  .  F  _  is  decreasing  in  n,  if  n,  >  2,  n0  >2, 

-  a.n^  ^2  lit 


a  <  .10. 
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Lemma  3.  Fi_a;n  >n  Is  Increasing  in  n2  if  n1  >  2 
n2  >  2,  a  <  .10  . 

Proof: 

Fa:nf,n2  <  Fa:nrn2  1f  nl  >  ni  (by  Lemma  2). 


So 


r1  <  f-1 

l-a:n2,n^  l~on2»n^ 


And  from  this  the  result  follows. 

% 

Using  the  Lemmas  in  expression  (3.24)  we  have 
1)  C2  <  Cj  since 


l-a:nr-  l-a^.ng  »-F«:n1  ,n2 


then 


-1 


-  1  ♦  A^,  +  BiF:'  -  1  +  A,F,  +  B,F 

C,  >  1  ,  C,  <  1  so  — U- - UL  < 

1  -  t  —  n,  +  c,f 


ri T  “ri1 


v  2 


nl  +  C2F2 


and 


so 


2>  f«!n,.  iFi]a:n,,  50  A1  i  A2 


-1  +  A]F1  +  B^1  -  1  +  A„F,  +  B^F;1 

n,  +  C,F,  < 


'2rl  Dlr  1 


'1  1  2 


n2  +  C2F2 


3)  F'1  F  >  F”^  F 

a:n^ ,«  a:n-j  .n^  -  1  -a : n-j  ,•  l-a:n^  ,n. 


Fa:n^ ,«Fa:n^ ,n^  —  "  Fl-o:n^ ,»Fl-a:n^ ,n^ 


o-F  f  „ )  <  o-f:1  F1  ) 

a:n-j .«  azn^.n^'  -  l-a:np« 


and 
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--1 


»«»Fa:nj  ,n,^Fa:nj »n^  —  ^  Fl-a:nj  ,®^l-a:n'j  *n^  l-am^.n^  V 


then 

and 


B1  —  B2 


-  1  +  A,F,  +  B^f;1  -  1  +  A2F1  +  B2f]] 


Yi  T1 
"l  +  C1F2 


nl  +  C2F2 


Going  back  to  (3.23)  we  get 

w  * '  -  (’-w)  ■  ('-w)  * ° 


w  -  w  *  w  - ' 


(3.25) 


Numerical  Integration  will  be  used  to  evaluate  the  integrals 

In  PL(0A)  »  PyU^).  anc*  P2^0A^’  method  use<J  simPson's  Rule 

which  gives 
% 

f^y)dy<^^  |fu)  +  ^f(a+h)  +  2f(a+2h) 

+  4f(a+3h)  +  2f(a+4h)  +  ...  +  2f(a+(n-2)h)  +  4f^a+(n-l)hj 

+  f(b)j 

where  h  Is  the  length  of  the  subintervals  In  which  (a,b)  has  been 
partitioned.  We  used  h  =  .001. 

Various  combinations  of  sample  sizes  were  cnosen,  and  for  each 
of  the  combinations  several  sets  of  ,  and  ,  and  different 
a's  were  used. 
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A  list  of  the  values  of  I,  J,  p^,  p^,  and  a  used  Is  given 

below 

(I.J)  •  (3»3);(3,7);(3,11);(5,3);{5,5);(5,11);(7,3); 

(7,7);(7,11);(9,3);(9,9);(9,11);(11,3);(11 ,7);(n,ll) 

oa?/o£2  :  .003,  1,  10,  100  (4.1) 

o  2/o  2  :  0,  1,  10,  100 
l  e 

1-a  :  .90,  .95,  .975,  .99 

V  1  +  J 

P2  "  1  +  I  OT2/°e2 

The  results  are  given  in  Tables  2,  3,  and  4.  The  first  five  columns 
in  each  table  show  the  I,  J,  n^ ,  n^,  and  n^  values.  Each  of  the  remaining 
columns  give  the  results  for  the  specified  1-a  confidence  coefficients. 
For  each  of  these  columns  the  ranges  that  the  confidence  coefficients  vary 
(as  p.j  and  p2  vary  over  their  allowable  values)  are  shown. 


The  approximate  confidence  coefficients  corresponding  to  the 
lower  confidence  intervals  are  very  close  to  the  specified  ones,  shown 
in  Table  2.  The  approximate  confidence  coefficients  for  the  upper 
confidence  intervals  range  from  values  very  close  to  the  specified  ones 
to  some  high  values,  shown  in  Table  3.  Similar  results  hold  for  the 
two-sided  confidence  limits  shown  in  Table  4. 

For  I  fixed  the  range  of  the  confidence  coefficients  is  closer  to 
the  nominal  value  when  J  increases.  For  J  fixed  the  upper  bound  of 
the  range  is  higher  when  I  increases.  So  the  confidence  coefficient 
will  be  better  when  J  is  large.  However,  in  general  the  approximate 
confidence  procedures  seem  adequate  if  I  and  J  are  each  5  or  bigger. 


Ranges  of  confidence  coefficients  for  lower  confidence  Intervals  on 
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Table  4 

Range  of  confidence  coefficients  for  a  two-sided  confidence 

interval  on 


I 

J 

nl 

n2 

n3 

1  - 

a 

=  .90 

1  - 

a 

=  .95 

3 

3 

2 

O 

k 

4 

.90103 

- 

.93208 

.95072 

_  _ 

.97000 

3 

7 

2 

6 

12 

.90043 

- 

.91349 

.95030 

- 

.96000 

3 

11 

2 

10 

20 

.89950 

- 

.90849 

.94827 

- 

.95702 

5 

3 

4 

2 

8 

.90096 

- 

.93896 

.95062 

- 

.97273 

5 

% 

5 

4 

4 

16 

.90077 

- 

.92684 

.95051 

- 

.96725 

5 

11 

4 

10 

40 

.89929 

- 

.91272 

.95031 

- 

.95884 

7 

3 

6 

2 

12 

.90094 

- 

.94310 

.95604 

- 

.97436 

7 

7 

6 

6 

36 

.89871 

- 

.92475 

.94984 

- 

.96607 

7 

11 

6 

10 

60 

.89896 

- 

.91673 

.94955 

- 

.96130 

9 

3 

8 

2 

16 

.90092 

- 

.94619 

.95059 

- 

.97647 

9 

9 

8 

8 

64 

.89835 

- 

.92374 

.94929 

- 

.96546 

9 

11 

8 

10 

80 

.89878 

- 

.92016 

.94938 

- 

.96334 

11 

3 

10 

2 

20 

.90090 

- 

.94863 

.95058 

- 

.97782 

11 

7 

10 

6 

60 

.89790 

- 

.93203 

.94926 

- 

.97006 

11 

11 

10 

10 

100 

.89868 

- 

.92315 

.94927 

- 

.96510 
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